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Abstract 

The  use  of  statistical  techniques  to  characterize  composite  materials  systems  and  their  use 
in  complex  electromagnetic  environments  is  discussed.  Examples  of  the  calculation  of  uncer- 
tainties and  distribution  functions  of  wave  statistics  are  given. 


1 Introduction 

Heterogeneous  media  which  consist  of  discrete  inclusions  of  various  kinds  and  configured  in  cer- 
tain arrangements  inside  a host  medium  can  be  considered  as  complex  electromagnetic  (EM) 
material  systems.  A rigorous  deterministic  analytical  or  computational  approach  to  their  anal- 
ysis is  of  almost  forbidding  complexity.  Hence  one  usually  takes  recourse  to  simplified  effective 
medium  theories.  Typically,  only  mean  values  (first-order  moments)  are  considered,  but  this 
tends  to  be  a significant  oversimplification  when  correlating  theoretically  predicted  characteris- 
tics of  a composite  with  measured  results  for  a realistic  sample.  In  particular,  variability  between 
different  realizations  of  samples  proves  to  be  a dominant  factor  in  the  uncertainty  budget  for  the 
macroscopic  constitutive  parameters.  The  resulting  uncertainties  are  furthermore  important  in 
assessing  the  accuracy  of  the  measurement  itself,  in  order  to  decide  on  the  significance  of  certain 
measured  effects  (e.g.  nonreciprocity  or  chirality),  or  on  the  relevant  truncation  point  in  any 
series  expansion  for  their  characterization  in  the  long-wavelength  regime.  In  crystal  physics,  the 
number  of  scattering  centers  is  usually  very  large  so  as  to  warrant  a negligible  level  of  uncer- 
tainty. However  for  synthetic  composites  containing  a relatively  small  number  of  inclusions,  this 
is  generally  not  the  case.  The  uncertainty  associated  with  sample  realizations  calls  for  new  or 
improved  methods  for  their  characterization.  Statistical  characterization  proves  to  be  an  espe- 
cially powerful  method,  because  the  uncertainties  decrease  with  increasing  degree  of  complexity 
of  the  material  system.  The  results  find  application,  for  example,  in  the  characterization  of 
adaptive  material  systems  in  active  or  passive  mode  of  operation  [1]. 

Secondly,  a growing  trend  exists  towards  the  use  of  complex  EM  environments  as  alternative 
measurement  techniques,  in  order  to  characterize  materials  in  their  operational  environment 
more  realistically  and  accurately.  Here,  realistic  excitation  and  illumination  conditions  are  be- 
ing generated,  as  opposed  to  the  idealized  case  of  single  plane-wave  illumination.  The  incident 
wave  must  then  be  considered  as  quasi-statistical,  for  the  measured  effect  is  due  to  an  ensemble 
average  of  a multitude  of  different  directions  of  incidence  or  polarization.  The  same  is  usually 
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true  for  the  internal  field  inside  the  medium  even  if  the  externally  incident  wave  is  deterministic. 
One  such  complex  EM  environment  simulator  is  the  NPL  stadium  reverberation  chamber,  in 
which  an  ensemble  of  illuminations  is  being  generated  for  the  EMC  testing  of  equipment  or 
characterization  of  EM  materials  or  systems  [2,  3].  The  most  general  case  is  of  course  the  one 
of  complex  media  subjected  to  complex  fields. 


2 Statistics  of  Waves  and  Media 


2.1  Moments 


In  the  general  case  of  complex  media  subject  to  complex  waves,  both  the  excitation  and  the 
constitutive  parameters  exhibit  statistical  fluctuations.  Sample  variability  can  be  taken  into 
account  by  incorporating  a continuous  or  discrete  realization  parameter  r.  For  example,  in  the 
Lorenz-Lorentz  formula  the  number  of  inclusions  and  their  dipolarizabilities  a thus  become 
random  variables.  At  a given  frequency  u> : 


e(r;u>)  = e0  + 


N(t ) a(r-,uj) 

1 - N(t)  a(r;w)/( 3e0) 


(1) 


The  average  of  the  macroscopic  electric  polarization  Pe  — Na/[  1 — Na/(3e0)]E  then  satisfies: 

(Na)E+N*(E)  + !?g-Pc 

V*e/  — , _ Na  vAJ 

L 3e„ 

where  (-)  signifies  ensemble  averaging.  Its  variance  then  follows  from  = (P2)  — (Pe)2 , 
which  also  takes  the  deviation  of  the  spatial  distribution  of  the  inclusions  from  a pure  random 
distribution  into  account.  With  D — eaE  + P*.  we  can  calculate  cr2D.  Its  expression  for  gen- 
eral bianisotropic  media  has  been  obtained  but  is  cumbersome.  For  the  case  of  an  anisotropic 
dielectric: 


= (£f  • (|t)  • <| )-E)  + (<Pt)-|t-|-<P))  +2  (#  • (|f)  • | • (E)')  - {E*-g t)  .(g.E)  (3) 
A similar  relation  exist  for  a\.  For  an  isotropic  medium,  (3)  can  be  simplified  to: 

al  = { D2)-(D )2  = a2eE2  + e2a2E  - 2eE(e)(E)  (4) 


where  ae  represents  the  standard  deviation  resulting  from  differences  in  material  processing 
conditions,  as  well  as  differences  or  uncertainties  in  the  number,  spatial  distribution,  size,  shape, 
etc.  of  the  inclusions. 


2.2  Distributions 

A full  characterization  of  the  wave  or  medium  ensemble  requires  knowledge  of  all  higher-order 
statistical  moments  or,  equivalently,  the  associated  probability  density  function  (pdf).  Provided 
the  pdf  of  the  constitutive  effective  medium  parameters  is  known  a priori  (which  is  often  the 
case  owing  to  the  central  limit  theorem),  wave  statistics  such  as  field  correlation,  wave  number, 
wave  impedance,  power  spectral  density  etc.  can  be  derived.  If  the  medium  can  be  considered 
as  ’random’,  which  presumes  that  the  number  of  inclusions  is  sufficiently  large,  the  ensemble 
permittivity  and  permeability  exhibit  an  approximately  Gauss  normal  distribution,  on  physical 
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grounds.  For  the  random  wavenumber  K or  random  refractive  index  N,  the  resulting  distribution 
is  then  obtained  upon  the  subsequent  variate  transformations  x — fie  and  n = \/x  as: 


/at(«)  = 


7T(TM0rc 


(5) 


where  K0(-)  is  the  modified  Bessel  functions  of  the  2nd  kind  of  order  zero.  For  the  wave  random 
impedance  Z,  the  transformations  y = fi/e,  z = ^/y  yield: 


fz(z ) = 


V~z 

* [M^e)z2  + (<TeM\ 


(6) 


For  the  mean-normalized,  statistically  isotropic,  homogeneous  and  unpolarized  field,  whose  three 
complex  components  exhibit  a circular  Gauss  normal  distribution,  the  real  (r)  and  imaginary 
(i)  parts  of  D and  B are  distributed  as: 


/^(cT1')  = (ncrecrE)  lKo 


(7) 


The  random  power  density  S associated  with  this  isotropic  field  satisfies,  for  deterministic  con- 
stitutive parameters,  a x%  pdf: 


fs(s) 


sp  1 exp(— s/2) 
2 p (p-  1)! 


(8) 


with  p = 1 for  a Cartesian  component  and  p = 3 for  the  total  rms  power.  For  the  field  magnitude, 
a X2 v Pdf  applies: 


f\Ea\(\ea\) 


\ea\2p  1 exp(— |ea[2/2) 
2 p-1  (p  — 1)! 


0) 


Other  wave  statistics  are  obtained  along  similar  lines,  starting  from  basic  normal  random  com- 
ponents with  a Gauss  normal  distribution. 


2.3  Spatial  dispersion 

Finite-size  effects  manifest  themselves  not  only  in  the  size  of  the  material  sample,  but  also  in 
that  of  the  inclusions  themselves.  For  random  fields,  the  effect  of  the  latter  on  the  effective 
properties  can  be  taken  into  account  by  applying  theoretical  results  for  general  local  averaging 
[2].  The  analysis  then  shows  that  finite  size  effects  give  rise  to  an  decrease  in  the  uncertainty  of 
the  constitutive  parameters,  but  to  an  increase  in  the  perceived  (i.e.  measured  ) randomness  of 
the  medium,  as  measured  by  its  normalized  spectral  bandwidth.  Typically,  for  sufficiently  small 
size  of  the  inclusions,  the  uncertainty  will  vary  according  the  square  of  the  characteristic  length 
of  the  inclusion,  or  with  the  square  of  the  averaging  distance  for  the  statistical  field. 
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